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Abstract 

We analyse the pairing symmetry of the order parameter in su- 
perconducting Sr2Ru04. Basing on a realistic three dimensional 3 
band energy spectrum we have introduced effective attractive electron- 
electron interactions and found a solution which has line nodes in the 
gap in a and (3 bands but not in 7 band. This state breaks time rever- 
sal symmetry and leads to the proper temperature dependence of the 
specific heat, in-plane penetration depth and thermal conductivity. 

1 Introduction 

Sr 2 RuO"4 is an oxide with a lattice structure similar to the well known high 
temperature superconductor La 2 _a;Ba x .Cu04 [l], with ruthenium replacing 
copper. When pure it becomes a superconductor at T C =1.5K and is one of 
the best candidates for spin triplet superconductivity || ||. [| • 

In spite of many convincing arguments for p-wave pairing the proper 
symmetry of three band superconducting order parameter is still unknown. 

The muon spin rotation experiments || indicate the existence time re- 
versal symmetry breaking, consistent with the order parameter A(k) = 
ia y a ■ d(k), with d(k) = (0,0,cf (k)) and d z (k) = A(T)(k x + iky). At the 
same time many thermodynamic properties show power low temperature de- 
pendences at low T, clearly pointing at the existence of nodes in the order 
parameter. The experiments on specific heat 0, penetration depth J?| and 
thermal conductivity || all show a quadratic dependence on temperature. 
This is surprising since a complete symmetry analysis of p-wave pairing in 
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tetragonal crystals shows no states which both break time reversal and pos- 
sess line nodes [|, . 

In the literature there exist a number of proposals for the relevant sym- 
metry states || H, [11], [12], [13], [14], [15|, [16]] . Most of them concentrate on two 



dimensional band models with p-wave symmetry states. Here we propose to 
start with the real three dimensional electron energy spectrum with all three 
bands (called a, (3 and 7) taken into account and allow for both in plane and 
out plane p-wave symmetries of the order parameter. 

In the next section we present our Hamiltonian and the approach, which 
is the full solution of the Bogolubov - de Gennes equations at each point of 
the three dimensional Brillouin zone. In principle our theory does not contain 
free parameters. By fitting the calculated Fermi surface to the experimentally 
determined one we fixed the single particle parameters of the Hamiltonian. 
The remaining three interaction parameters have been fitted by requiring 
that all three bands have a single superconducting transition temperature 
(i.e. the monotonic low T behaviour of the specific heat) and a correct value 
of the transition temperature. Using the parameters we have solved the 
Bogolubov-de Gennes equations. The obtained gap functions have been used 
to calculate penetration depth and thermal conductivity. 



2 The Symmetry of the Order Parameter 

We start with the description of the electronic structure of Sr 2 Ru0 4 in the 
vicinity of the Fermi energy ep in the normal state. We consider three orbitals 
of ruthenium A, B, C of character d xz , d yz and d xy , respectively. These 
generate 3 Fermi surface sheets a, (3 and 7. The a and /3 sheets arise from 
hybridised A and B orbitals which we model with hopping parameters (t±, £2, 
Iab), which describe AA (BB) hopping in x(y) direction, AA (BB) hopping 
in y(x) direction and AB hopping, respectively. The 7 sheet arises form the 
C orbitals, which we model by hopping parameters t and t' for nearest and 
next nearest neighbour sites. We have fitted these parameters to the de Haas 
- van Alphen experimental cross-sections of the Fermi surface flT7 |. To allow 



for the three dimensional Fermi surface |18| we add the hopping parameter 
t± which we assume to take on the same value for each band. 

Due to the symmetry of the orbitals and the body centered tetragonal 
symmetry of the crystal we assumed that the dominant pairing interactions 
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Figure 1: Eigenvalues E l3 '' y (k) at the (3 (a) and 7 (c) Fermi surfaces plotted 
as a function of k z and xy plane azimuthal angle. 

were in plane C/ 7 for the 7 band, and inter-plane C/ Q)( g for the a and /3 bands. 
Each of these was tuned by fitting to the experimental superconducting crit- 
ical temperature Tc = 1.5 K, so we have no free parameters. 

Our strategy is to solve the Bogolubov-de Gennes equations to find the 
pairing potentials and energy spectrum. We allowed for the following general 
structure of the pairing potential 



+ 
+ 



A^'sin^ + A^'sin^ 

A u x f sin (k x /2) cos (k y /2) + A v y f sin (k y /2) cos (k x /2)] cos (k z c/2) 

(1) 



A u z f sin (k z c/2) cos (k x /2) cos (k y /2), 



1, and where indices u, v! take the 



in units where the lattice constant, a 
values a, (3, 7. 

For the given set of parameters the only nonzero pairing potentials have 
the following symmetry AT = iA y ; a , A^f = iA 13 / A™ = iA™. Each of 
these solutions therefore breaks the time-reversal symmetry. Moreover, due 
to the factor cos(k z c/2) (Eq. 1) we have planes of zeros on the a and (3 sheets 
for k z = ±71 /c. 

To illustrate the structure of the gap on the Fermi surface we have plotted 
in Fig. 1 the eigenvalues for the (3 and 7 sheets. Figure (la) illustrates clearly 
the zero eigenvalues in the f3 band on the planes k z = ±7r/c. The structure 
of the gap on the a sheet is very similar. Contrary to that the gap on the 7 
band (Fig. lb) is always finite but has strong modulation in the (x,y) plane. 
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Figure 2: (a) Pairing potentials |A| 

bands a, (3, 7 as a function of temperature, (b) The specific heat C/T versus 
temperature (the solid line corresponds to the calculated results, while the 
points are the experimental data H). 



The temperature dependence of the gap parameters is shown in Fig. 2a 
and the specific heat in Fig. 2b. The overall agreement with the experimen- 
tal data is good. In particular the height of the step at T c and the linear 
dependence of C/T on temperature agree with experiment 0. 

For the above state we have also calculated the penetration depth A (Fig. 
3a) and the thermal conductivity k (Fig. 3b). The penetration depth shows a 
close to linear dependence on temperature at low T and agrees approximately 
with the experimental data . Note that impurity scattering or other effects 
may change the low temperature power low. The thermal conductivity is 
quadratic as expected for a gap with line nodes. 

3 Remarks and Conclusions 

Using a realistic three band model of S^RuC^ we have found a p-wave 
pairing state which has nodes and breaks time reversal symmetry. The state 
has 3 -E'ji[e] (axial) pairing symmetry in the notation of |§. The line node on 
the a, (3 sheets is a consequence of the specific form assumed for the pairing 
interaction. 

The state we have found is of similar nature as the one proposed on 
different grounds in [15|, |l(J. Our model contains only attractive interactions 
unlike one proposed in ||15|| . The nice feature of our model is that we have no 
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Figure 3: The in plane penetration depth A for low temperature (a) (the 
diamonds correspond to the results of calculations while the solid line is a 
linear fit to them). The thermal conductivity k as function of temperature 
(b). 

adjustable parameters. The state we obtained leads to correct temperature 
dependence of not only specific heat but also the penetration depth and 
thermal conductivity. 
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